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COMPARISON OF FINITE SOURCE AND PLANE WAVE 
SCATTERING FROM CORRUGATED SURFACES 

D. M, Le Vine 
ABSTRACT 

The choice of a plane wave to represent incident radiation in 
the analysis of scatter from corrugated surfaces is examined. 
This is done by comparing the physical optics solution ob- 
tained for the scattered fields due to an incident plane wave 
with the solution obtained when the incident radiation is pro- 
duced by a source of finite size and finite distance from the 
surface. The two solutions are equivalent if the observer is 
in the far field of the scatterer and the distance from observer 
to scatterer is large compared to the radius of curvature at 
the scatter points, conditions not easily satisfied with ex- 
tended scatterers such as rough surfaces. In general, the 
two solutions have essential differences such as in the loca- 
tion of the scatter points and the dependence of the scattered 
fields on the surface properties. The implication of these 
differences to the definition of a meaningful radar cross sec- 
tion is examined. It is shown that the radar cross section 
defined from incident plane waves is meaningful for the case 
of finite sources and extended scatterers if the far field 
conditions are met or, if the scatter is incoherent and the 
surface statistically homogeneous, whenever the observer 
is far from the surface compared to its radius of curvature 
at the scatter points. 
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COMPARISON OF HNITE SOURCE AND PLANE WAVE 


SCATTERING FROM CORRUGATED SURFACES 

The choice of a plane wave to represent incident radiation is a common and use- 
ful approximation in scattering problems. On the other hand there are many 
situations in active (radar) sensing of the earth from airborne and space plat- 
forms where the specific phase structure of the incident radiation can be im- 
portant, Also, in cases of extended scatterers there is no a priori guarantee 
that radar cross section as derived from analysis of plane waves is applicable 
for predicting scattered power. It is the objective of this article to compare the 
scattered fields obtained by assuming an incident plane wave with those obtained 
by assuming a finite source for the case of scatter from stochastic surfaces. 
Some obvious differences are manifest, such as amplitude decrease due to wave 
spreading, present in the case of radiation from a finite source but not in the 
case of plane waves, and differences in the temporal history of the received 
signal which are a consequence of the plane phase structure of the incident radi- 
ation with plane waves and the curved structure with finite sources. More 
subtle differences involving the dependence of the solutions on the parameters 
of the surface also appear, and may have consequences in interpretation of 
scattering cross section. 

The specific example to be considered here is scattering from irregular (i. e, , 
stochastic) surfaces with the restriction to perfectly conducting boundaries 
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and two dimensions. This is an idealized model relevant in first order for 
scattering from ocean surfaces and relevant to some extent for scattering 
from plowed fields and perhaps rough snow and ice. The restriction to two 
dimensions has been chosen for the obvious simplification in treating the 
vector problem. A physical optics solution (Kirchhoff approximation plus and 
asymptotic evaluation of the Helmholtz integral) is obtained, and is a reason- 
able approximation for microwave radar scattering near nadir, at least for 
ocean surfaces (BasSj et al. , 1968; Barrick, 1968a-b). The object of the 
analysis is to compare the fields scattered back to an arbitrarily located 

observer when, in the first case, the incident signal is a plane wave, e(r, u) ~ 
ik* 7 

E^e'' , and in the second case in which the finite nature of the radar antenna 

is taken into account. The geometry is illustrated in Figure 1. 

In the following sections the scattered fields will be computed for the two cases 
and a comparison made. Then, a radar cross section applicable to extended 
scatterers will be defined from the radar equation, and the solutions for the 
scattered fields will be used to identify conditions under which this definition 
of cross section corresponds to the conventional definition obtained from 
incident plane waves. 
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SCATTERED FIELDS 


Both solutions for the scattered electric field, e_(r, v), can be written as fol- 

s 


lows in the case of perfectly conducting boundaries: 




\v)= - ^ g(r/r')ds 

surface 


( 1 ) 


where g(r/F’) is the two-dimensional Green’s function, j/4 -71), 

and n is a unit vector normal to the surface, Z{y). Assuming perpendicular 
polarization {E^ = ^ ^ Kirchhoff approximation yields the following 

result in the case of plane waves: 


3 e 
dn 


2 j (k- n) T(i') Eq X 



( 2 ) 


where j(i^) is the Fourier transform (temporal) of the incident pulse. In the 
case of a finite source, assuming that the source can be represented in 
terms of an equivalent current distribution, x j(i') J(x, y), the Kirchhoff 
approximation yields the following result (Appendix A) : 


source 

where R = yjiy-y')^ ^ t} (y’)* find (kR) is the derivative with respect to 

kR of the Hankel function. 


j(i') 


/ 


k^z(y') ( 1 )' _ 

J(y'. z') H; ' (kR)dr 


2 R 


( 3 : 
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The scattered fields are obtained by substituting Equation 2 or 3, as appro- 
priate, into Equation 1. Performing the Integration in the limit that kR— 
by means of a saddle point approximation (Copson, 1971), yields the following 
result for backscatter in the case of plane waves: 


ep(r, t') = 


y cosie-c) jk*(y„). 

° 4- C08(e + a) 


all y 


2 COS 6 


1 - 


R(y„) 

COS Q cos^(6 + a) ^c^^n^ 


(4) 


where 6 = tan"* 


kz 


ky 


-1 r^zi 

is the angle of incidence; a(yn) = tan 
is the slope of the surface; R(yn) is the distance between the scatter point 
at and the observer; R^iy^) i® the radius of curvature of the surface at the 
scatter point, and 


<t>(yn) = yn "Z(yn) cos 6 + R(yn) 

The scatter points, y^, are determined by the requirement that 
and for backscatter the requirement 0 < o pertains. 

In order to obtain the equivalent expression in the case of a finite source. 
Equation 3 must first be evaluated. For this purpose, assume that the surface 
is far enough from the source to permit the fraunhofer approximation to be 


(5) 


94> 

Ty" 


= 0; 
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made in Equation 3. Then, using an asymptotic form for (kR) and using 
a binomial expansion for R one obtains (Appendix B) : 


be 

dn 


1/2 V^/t' j(»/)k^ cos(v>-a) (kR) F (y, ) 


( 6 ) 


where R is the distance from the center of the source to the surface and 

F (y,v ), is the Fourier transform of J(yz) evaluated at spatial frequencies, 

M = ~ sin <P and cos ^ where 

y c z c 


(P (y) = cos"‘ 


H - Z(y) 
R 


(7) 


H is the height of the source above the mean surface (i. e. , the ordinate of 

( 1 )' 

the source). The tilda (~) above (kR) in Equation 6 denotes the asymp- 

totic form for large kR, and the prime indicates the derivative with respect 
to R: (kR) ^ jk Vz/kR exp [j (kR - ff/4)j. Substituting Equation 6 

into Equation 1, and evaluating the integral in the limit of large kR, one 
obtains: 


ejf, i>) = X n/pA 


j(i^) V k/16 



ei2kR(yn) . 




(B) 


Equation 8 is the electric field scattered back to the finite source. 
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In order to compare Equations 4 and 8, it is first necessary to find an equiv- 
alent in the case of the finite source for the amplitude of the plane wave. For 
this purpose assume that the amplitude. £^. of the plane wave, is the magni- 
tude in the fraunhofer limit of a cylindrical wave radiated by the finite source. 
Since the far field radiated by the finite source in the direction of incidence 
of the plane wave (i. e, , 0 ) is : 

e(f, v) = jkcxAj^ 

rr~, — l f . jk [y sin o -zcos 0 ] , , 

= -x\M/e j(v)\k/87r' e ^ * I J(y,z)e L-" Jdydz 

J 

- oe 

^jkR 

= F(0.v) (9) 


one obtains: 


eJ(m) = - VmA ' e •’ Vk/8ff j(v) F(0, v) (10) 

Using this definition, one obtains the following formulas for the scattered 
fields: 


a) Plane wave 


e (r.v) = - EJtv) x 




cos (e?-g) jk4>(yji) 
cos (0 + Q ) 


2 cos e 


1- 


R<yn) 

COSCr co8*(0+a) R-ty^) 


1-^ 


(lla) 
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b) Finite source 


eg(r,t>) = - gj(i^) 

ally. 




1 - 




(lib) 


where F{y^,^)= F(y^j, t- )/F(0, i-). If F( 0, v ) is the maximum value of 
F(yni v) then Fiy^^, v ) Is the relative field pattern for the source {Collin and 
Zucker, 1969). Also the following definitions have been made assuming the 
observer to be at (0, H): 


4>(yn) = R(yn^ yn ^ “ ^^yn^ ® 

(12a) 

R(yn> = V[H-Z(yj^) + y/ 

(I2b) 

EqTi v) = - e ^ k/8ff' j( v) F( 0, v) 

(12c) 
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i OMPARISON 


In comparing Equations 11a and 11b it is important to keep in mind tha t 
Equation 11a represents the scattered fields when an ijicident plane wave 
(source at infinity) scatters to an observer at some finite distance from the 
surface; In contrast. Equation 11b represents the fields in the case of back- 
scatter from a source of finite size and at finite distance from the surface. 

In the former case, plane phase fronts scatter (in the physical optics limit) as 
cylindrical waves back to the source whereas in the latter case, incident 
cylindrical waves scatter as cylindrical waves back to the source. This dif- 
ference accounts for the additional distance factor, \ R^^, in the denominator 
in Equation lib; it represents cylindical spreading of waves on transmission 
from the finite source. 

Other amplitude differences as well as differences in phase and location of 
the scattering points also occur. For example, consider the location of the 
scatter points, y^^. The physical optics solution requires that in both cases 
these be "specular” points — angle of incidence equals angle of reflection 
relative to the local surface normal — but the geometrical possibilities are 
much different in the two cases. In the case of an incident plane wave all 
incident rays have the same direction whereas for the finite source they are 
all radial vectors emanating from the source. As a consequence, in the case 
of the finite source, backscatter requires incident and reflected rays to 
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coincide and to be perpendicular to the surface at the scatter point. In con- 
trast with au incident plane wave all manner of variations are possible includ- 
ing what should rightly be called forward scatter. (See Figure 1. ) As a re- 
sult the distances, and R(Yj^), are not the same in the two solutions. 

Only if one restricts the size of the illuminated surface will 4 > (y^^) , R(y^) and 
the location of the scatter points, y^^, be compa 'ble in the two cases. For 
example, if one were to require that the observer be in the ’’far field" of the 
illuminated surface, then expanding R(y) in a binomial series about the 
distance R^ from observer to center of the illuminated footprint, one obtains: 

R(y) + y sin 0 - Z(y) cos 9 (13a) 

and consequently 

4»(y) ^ R(y) + y sin 0 - z(y) cos B (13b) 

^ 2R(y) - Rq 

That is, except for the constant phase factor, kR^, which is arbitrary for the 
plane wave, the two phase factors are essentially equal when the scattering 
surface is small. Consequently, in this case the location of the scatter points 
will be the same for either type of incident variations. 

But even with the far field restriction, the amplitude terms in Equations 11a 
and 11b are different. That is, even with 0 = a in Equation 11a and assuming 
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the scatter points are the same, there is a fundamental difference in the 
radicals. This is again a manifestation of the different phase structure — plane 
versus cylindrical — incident at the scatter point in the two cases. Only if th^ 
distance from surface to observer is much greater than the radius of curvature 
at the scatter point (i. e. , R(yn) » Rc(yn)) radicals become equal. This 

is because in the physical optics solution the incident signal scatters from an 
equivalent mirror at the surface with focal length, Rc(yn)» and the phase dis- 
tribution along this infinitesimal mirror is functionally different in the two 
cases treated here. 
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IMPLICATIONS FOR RADAR CROSS SECTION 


The restrictions mentioned above are sufficient to insure that Equations 11a 
and 11b are essentially term-wise equivalent (with appropriate consideration 
given to antenna gain and range effects). That is, if the observer is in the 
far field of the scatterer and sufficiently removed that the distance to the 
scatterer is much greater than the radius of curvature at the scattering points, 
then Equations 11a and 11b differ only in factors due to the radiation pattern 
of the antenna and cylindrical spreading of the incident radiation which are 
present in the case of the finite source. 

The interesting case occurs for extended scatterers such as irregular sur- 
faces where one can be far from the surface but receive energy scattered 
from a wide range of incidence angles. Such situations arise, for example, 
in monitoring ocean surfaces either with an altimeter where the pulse first in- 
tersects the surface at nadir and then at increasing incidence angles as a func- 
tion of time (Bar rick, 1972), or in such off nadir sensing configurations as is 
employed with the Short Pulse Radar where a long narrow antenna beam 
directs the pulse along the surface from an initial intersection near nadir 
to an eventual intersection at 30° or more (Tomiyasu, 1971; Le Vine et al. , 
1975). 

In case of extended scatterers one might imagine computing the scattered power 
as a sum over the power scattered by rays incident at the various angles allowed 
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by the scattering geometry. For example, based on the radar equation, one 
might write (in two-dimensions) the following form for the received power; 

r 4 Pt(v)Gt(0)GR(0) 

J T 57^5 (14) 

surface 

where a°{o} is a cross section per unit surface (length in this two-dimensional 
example). A three dimensional equivalent to this equation has been used to 
analyse the performance of a radar altimeter (e. g. , Harger, 1972) and has 
been used traditionally in interpretation of radar scatter from the ocean (Kerr, 
1951). (See Appendix C for definition of terms. ) 

o 

Equation 14 can be regarded as a definition of a {&): That is, one determines 
0 ° (0) by formally computing the received power using solutions for the scat- 
tered fields (such as Equation 11b) and putting the results in the form of 
Equation 14. In two important cases this procedure leads to the same cross 
section as obtained with the definition bs sed on incident plane waves (Kerr, 
1951; Skolnik, 1970). These cases are: 1) When the observer is in the far 
field of the scatterer and far away compared to the radius of curvature at the 
scatter points; and 2) When the scatter is incoherent, the surface is statis- 
tically homogeneous (i. e. , spatially stationary) and the observer is far from 
the surface compared to the radius of curvature at the scatter points. 
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The equivalence in the first case is a consequence of the discussion in the 
previous section. In this case the scatterer can be regarded as essentially 
a single point and the conditions of the previous section are adequate to insure 
that the two solutions, Equations 11a and 11b are identical. 


In order to demonstrate the equivalence in the second case expressions are 
necessary for the radar cross section from Equation 14 and as defined by the 
assumption of incident plane waves. The expression for the radar cross 
section as obtained from plane waves is obtained by a straight forward sub- 
stitution of the plane wave solution, Equation 11a, into the definition of cross 
section (Kerr, 1951; Skolnik, 1970). With minor generalizations to specialize 
the results to two dimensions, pulses and stochastic scatters, one obtains: 


( o°(d)) =— lim 2 7rp 
P ^ P-^oo 


e„(r, t/)i e„(r,v) 


EqJ( v)*E^j(t>) 


> 


(15) 


where the brackets ( > denote an ensemble (statistical) average and the 
asterisk-dot ( i ) denotes a convolution of the pairs in a scalar (dot) product. 
(The convolution is the frequency domain equivalent of a product in the time 
domain and reduces to a simple product in the case of harmonic time depend- 
ence and a complex analytic representation of the fields. ) Lis the length of 

the illuminated surface. The scattered field, e^^(r, i;), is given by Equation 

s 

11a. In the special case R(y ) » R (y ) and incoherent scatterers. Equation 

ii n 

15 becomes: 


14 


(16) 


<«;(«)> 



w 


gjk[v^(yn) - 



> 


where is an arbitrary reference phase. When the surface is homogeneous, 
this solution reduces to; 

where t? is the number of scatterers per unit length. 


The equivalent expression for distributed scatterers is obtained by computing 
scattered power from Equation lib and putting the result in to the form of 
Equation 14. The received power is obtained from Equation 11a as follows: 

G„ (r, „) 

where is the equivalent "area" of the receiving antenna and G„(r, v) 

k/4 

its gain(Friis and Lewis, 1947). Assuming that R{y^)» Rgiyn) and the scatter 
is incoherent, one obtains (Appendix C) : 





4 Pt('')GR(yn)Gt<yn)r 

^ l2RR(y„)l^ [”'W 


[R(yn)-Ro] 



Assuming homogeneous surface statistics and neglecting the small dependence 
of the antenna parameters (and distance) on the surface heights, this 
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expression can be written: 




f 4 Pt(>') r [R(yn)-Roll 

° J T (2,r)4 ^ [-4< R,(y„)e ^ 


surface 


( 20 ) 


Comparison with Equation 14 yields the following form for <«°(0)> : 


, I 1 j2k rR(yn)-Rol 

(a (0)> = rrn <|Rc(yn)|e J 


> 


( 21 ) 


which is the same as Equation 17 and is the two dimensional equivalent of re- 
sults which have appeared in the literature on scattering from ocean surfaces 
(Barrick, 1968; Kodis, 1966). 


The Important point to be made here is not so much the equivalence of Equa- 
tion 17 and 21, but rather that the equivalence is the consequence of several 
assumptions (far field, R(y^) incoherence and homogeneity). These 

assumptions are, of course, not universally true. It was shown above that 
the expression for the scattered fields, assuming incident plane waves and 
waves incident from a finite source, are not identical unless the observer is 
in the far field of the surface and much further away then than any significant 
radius of curvature on the surface. It follows, that the plane wave expression 
for cross section, o°(fl), is not always applicable for use in Equation 14. 
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CONCLUSIONS 


It has been shown that in order for the scattered fields obtained using a 
finite source and incident plane waves to agree, two conditions must be sat- 
isfied; 1) The scatter must be of finite extent with the observer in the far 
field of the scatterer; and 2) The distance from scatter points to the observer 
must be much greater than the radius of curvature of the surface at the scatter 
point. These conditions guarantee that the solutions differ only in the obvious 
spreading due to the cylindrical nature of the incident energy in the case of a 
finite source. If the preceding two conditions are met, and one properly 
keeps track of the spreading, then a radar cross section can unambiguously 
be defined using the definition based on plane waves. Otherwise, the equiv- 
alence of radar cross section as defined for incident plane waves and from 
the radar equation only pertains in special cases. One such special case occurs 
for incoherent scatter from statistically homogeneous surfaces if the observer 
is far from the surface as compared to the radius of curvature of the surface at 
the scatter points. 
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APPENDIX A 


FIELDS DUE TO A SOURCE ABOVE 
A PERFECTLY CONDUCTING PLANE 

Consider the two dimensional problem of finding the electric field due to a 
current source which is above a perfectly conducting plane and directed par- 
allel to the surface. In the case of an x-directed current, the x-component 
of magnetic potential Is adequate to describe the fields, and the electromag- 
netic fields are given by: 


e(r, »^) = 

jckA^lr, i') X 

(Al) 

h(F, v) = 

jF X [Ax(F,i^) xJ 

(A2) 


The magnetic vector potential Aj^(r, i>) can be expressed in terms of the two 
dimensional Green’s function, gj^(F/r'), as follows: 

Ax(r, K ) = J Jjj(r,' V ) g^(r/r ' ) dr ' (A3) 

The Green’s function can be represented as a sum of a free space (no bound- 
aries) Green’s function plus a scatter term due to the boundaries: 

gj^(r/r ') = g^(r/r') + g® ”') ^^4) 

and the scattered term can be obtained by finding the plane wave representa- 
tion of g® (F/r ') and then constructing g®(F/F') out of the same superposition 
but using the scattered (reflected) plane waves (Clemmow, 1956). Doing so. 
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one obtains: 


gj(r/r') = ~ Hj”(k7) 
g®(r/r ') = - ■“ Hi'’(k7 ) 


(A5) 

(A6) 


where for a source at (y'; z') (Figure 2) one has: 

7 = V(y - y')2 + (z - 

> = \/17- y')^ + (z + z')2 

The scatter component is just the free apace Green's function due to a point 
source at z' units below the surface with current in the negative x-direction. 
Using Equations A5 and A6, one obtains the following result for the x-compo- 
nent of the magnetic vector potential: 

= T f (kn - (k'T)] dr' (A7) 


Substituting this expression for A^(r, v) into Equation Al, one obtains 


Cjjtr.p) = - k/4j Jj^(r>) (k7) - (k7)] dP (A8) 

It is trivial to check that this expression is zero in on the surface (i. e. , at 
7 = T). The derivative of Cjjtf, u ) normal to the surface is also easily com- 


puted: 

dejj(r, V ) 

3z 


f 


z - z' 


z + z' 


kV4 j Jjj(r> 


Ho”’ (k7)--:r- Hq'*' (k=T) 

7 7 


(1) 


dr' (A9) 
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and when evaluated on the surface yields: 


de^ir.v) 

dz 

z 


Figure 2. 


0 




(A 10) 



Geometry for the Calculation of the Fields Due to a 
Source Above a Perfectly Conducting Plane 
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APPENDIX B 


INTEGRATION OVER THE SOURCE COORDINATES 


It was shown in Appendix A that in the case of a line source, ^ove 

a perfectly conducting plane, the derivative of electric field intensity normal 
to the surface and evaluated on the surface (z - 0) is: 


3ex(r, v) 


J (r',i^)Ho‘‘^ (k» df' 

27(f') 



The coordinates in Equation B1 refer to the configuration shown in Figure 2; 
however, because this analysis is to be applied to a plane tangent to the 
stochastic surface, Z(y), at many different points, it is necessary to express 
all coordinates in terms of a reference coordinate system (y, z) as shown in 
Figure 3. The coordinates in which Equation B1 is presently expressed are 
the primed system in Figure 3. The origin of the primed system has been 
chosen so that the z' axis passes through the center of mass of the source 
because it is convenient to do the integration over the source coordinates in a 
center of mass system. Letting the center of mass system be (tj, { ) with axes 
parallel to those of the reference system. Equation B1 becomes: 
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P'igure 3. Ixical Coordinate System for Calculating the Fields on the Surface (Kirchhoff Approximation) 



(B2) 




(Z'(^,T?)] 


(k7) dTJ dt 


where: 


7(17, U = ^fp(y) cos <iP + + [P(y) sin 

Z\v,^ ) = P(y) cos (<5 - a ) + ( ^cos a - t? sinor J 
p (y) = ^[R cos e - Z{y)P + { R sin e + y]2 


0 (y) = tan 


-1 


[ R cos © - 
R sin 8 + 


- ZJy) 

y 


Although this integration can not be done in general, it can be approximated 
in the important special case in which the source is in the far field of the sur- 
face in terms of the distance, P (y), between its center of mass and the point 
at which the plane is tangent to z(y). In this case, employing an asymptotic 
form for the Hankel function. Equation B2 oecome,' ; 


9e. 


X 


dz 


^ IZ 

ye 2 


(kp) I /cos (<P-a). 


■ ff. 


(B3) 


lz = 0 


JJ 




j 2 Jrfa(y)ti+ b(y)^j 


dr) d^ 


where 

V 

a(y) = - — sin [^(y)] 
b{y) = cos [>P(y)l 
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The limits integration may be formally extended to include the entire n , t 
space, and then one obtains: 

= (kp) cos (<P -oi) G(y. v ) (B4) 

2 

z = 0 

where G(y, v ) is a Fourier transform of the spatial coordinates of the source, 
J( I evaluated at frequencies tfjj- a(y) and = b(y). Whenever the 
spatial and temporal dependence in the source are separable [i. e. , J(r, if) = 
Jp(r) j(i')] one obtains: 

= j(i^) fj [Jr(T7,n] p^ = a(y) (^5) 

= b(y) 

where denotes a Fourier transform. 

An important special case is that in which ) is a delta function (i. e, , 

one has a point source). In this case the integrations are easily evaluated to 
yield: 

G(y, «/) = )(»») (B6) 

It is convenient in employing these results to separate explicitly the spectrum 
of the current waveform, and the spatial character of the source. Thus, 

let F(y, ) be defined by: 

G(y.t') = j(i') F(y, i') (B7) 

It is the factor, F(y, v), which is employed in the text. 
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The gain, of the transmitting antenna is defined to be the power radiated 

per unit surface in the direction, divided by the total power radiated per 
unit surface (ColiLi and Zucker, 1969). That is; 




1 

^ f F2(«P)d«^ 


(C3) 


The power received by an antenna is equal to the incident power density 
weighted by the effective area of the receiving antenna (which constitutes a 
definition of "effective area"). This effective area can be related to antenna 
gain by means of the reciprocity theorem. In the case of two dimensions the 
effective area of an antenna with gain, Gj.(r, p), is (Friis and Lewis, 1947): 


A^(r, y) = (4A)Gj.(r, v) 


(C4) 


The received power is given in terms of the effective area, Aj.(r, t^), and the 
scattered electric field, e„(r, i^), by: 

eg(r, ) \jA^(r,u) 


6g(r, i') \lA^(r,v) 


(C5) 


The form required in the text for received power scattered from the surface is 
obtained by substituting Equation 11b for the scattered fields into Equation C5 
above. Using the preceding definitions Equation lib takes the form: 


6g{r, i^) = X ^ 






yj2R{y^) 


2TT 

Gt(<P) f 

1 f 2(.^) d<P 


2tt 


R(yn) 


1- 


L ^c^^nU 


'j2k [R(yn)'-Rc(yn)] 


(C6) 
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where 


Iq(v) = i(v)yjl:/6n exp 


j [ 2k - V4 ] 


(C7) 


R (y ) is the distance from the observer to the mean surface at the scatter 
on 

point, y^. It has been introduced here to reference the time of reception 
of the scattered pulse to the round trip propagation time between source and 
surface. In this way, fluctuations about this mean time due to surface ir- 
regularities are included explicitly in the factor: exp j2k [R(7j^) - ^ 0 ( 7 ^) ^ 
Now assuming that R(yj^)> >Rj,(yj^) and that the scatter is incoherent, one 
obtains the following expression for the mean received power: 

21T 


4/k<[yv)::--Io(v) r F2(<P)d^] 

>*11 J. 


(2ttR)2 
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Rc(yn)l e 


j2k [R(y^) - RgO'n)] 


_ 




all y. 




(2ttR)^ 


j 2 k [R(y„) -Ro(yn)] 

Tr|Rc{yn)l e -* > 


(C8) 


which is the form used in the text. 
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